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Abstract 

We analyze the nature of the rotational energy level structure of the OCS-He and 
OCS-H2 complexes with a comparison of exact calculations to several different dy- 
namical approximations. We compare with the clamped coordinate quasiadiabatic 
approximation that introduces an effective potential for each asymmetric rotor level, 
with an effective rotation Hamiltonian constructed from ground state averages of 
the inverse of the inertial matrix, and investigate the usefulness of the Eckart condi- 
tion to decouple rotations and vibrations of these weakly bound complexes between 
linear OCS and ^He or H2. Comparison with exact results allows an assessment of 
the accuracies of the different approximate methods and indicates which approaches 
are suitable for larger clusters of OCS with ^He and with H2. We find the OCS- 
H2 complex is considerably more rigid than the OCS-He complex, suggesting that 
semi-rigid models are useful for analysis of larger clusters of H2 with OCS. 
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1 Introduction 



Until recently, Helium-4 has been the only known boson superfluid liquid. 
Many aspects of superfluidity have been studied experimentally, phenomeno- 
logically, and theoretically. These include its relation to Bose-Einstein conden- 
sation, the temperature dependence of the superfluid fraction in the 2-fiuid 
model, the finite size-dependence and the influence of atomic and molecular 
impurities in ^He on the superfluid fraction. The latter two issues have raised 
new, interesting and more basic questions. In particular, what is the validity 
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of phcnomcnological models on a microscopic scale, and how should we de- 
tect superfluidity of small systems. These questions ultimately demand a full 
microscopic understanding of superfluidity. 

Rovibrational IR spectroscopy experiments with chromophores in helium 
droplets [1,2] have conflrmed superfluid behavior for systems as small as 
approximately 60 '^He atoms [3] by observing the free rotation of the OCS 
molecule. Subsequently, a new flcld of research has been opened by similar 
experiments in small para-hydrogen clusters embedded in low temperature 
hehum droplets. [4] Spectroscopic experiments on OCS in these hydrogen clus- 
ters show spectroscopic anomalies which have been interpreted as indicating 
the existence of superfluidity of para-H2 for clusters small enough to be fluid- 
like rather than forming a solid shell around the chromophore. These conclu- 
sions have recently been verified by path integral calculations showing that an 
anistropic superfiuid fraction appears at low temperatures (T < 0.3 K).[5] 
In addition, new experimental studies on very small complexes with less 
than one solvation shell of '^He or H2 obtained rotational/rovibrational ex- 
citation energies for complexes of well-defined size, 0CS-pH2 and OCS-^HcAr, 
N — 1,. . . , 8.(6,7,8,9,10] New theoretical work for small OCS-^Hejv complexes 
[11] show an interesting transition from semi-rigid rotation of the whole clus- 
ter for small N to partial decoupling of the ^He motion and the OCS rotation 
for increasing A^. As explained in detail in Ref. [11], this implies a transition 
from a molecular complex to quantum solvation. 

Quantitative microscopic calculation of the spectrum of excitation energies of 
0CS-(pH2)Ar and OCS-'^Hcat clusters is an important step towards the under- 
standing of the rovibrational spectra of these unusual clusters. For ^ 1, 
this is still an extremely demanding task with current methodological and 
computational capabilities. Exact excitation energies are often not accessible 
and it is therefore necessary to make use of dynamical approximations, based 
whenever possible on existing knowledge of e.g. small cluster structures or 
on microscopic understanding of quantum solvation structures. In this paper, 
we therefore present a survey of approximate methods for the calculation of 
rotational energies of small Van der Waal complexes that have potential ap- 
phcation to large quantum clusters of OCS in ^He and in H2. In order to be 
able to critically assess the various approximations by comparing with exact 
results, we focus here on the smallest case, N — 1. 

The methods we employ here are all based on diffusion Monte Carlo (DMC) 
sampling. Our focus on DMC methods derives from their suitability for ap- 
plication to larger sizes that are beyond the scope of basis set approaches. 
DMC yields exact results for ground state energies and the approximations 
employed here all take advantage of this. The approximate methods all assume 
some degree of rigidity, an assumption which we expect would progressively 
fail to be justified as one would go to larger clusters, with helium failing at 
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smaller cluster size than hydrogen. In section 3 we give an overview of how to 
calculate approximate rotational excitations. We described the clamped coor- 
dinate quasiadiabatic approximation of Quack and Suhm [12], a modification 
of it which incorporates the Eckart condition [13] for semi-rigid complexes, 
and a perturbative calculation of rotational constants following Ref. [14]. In 
section 4, we compare these results with the exact bound state energies deter- 
mined by the program BOUND [15]. We discuss here also the relative rigidity 
of the OCS-He and OCS-H2 systems as apparent from the grouping of exact 
energy levels into (ro)vibrational and rotational subsets, where "vibrational" 
from now refers to the radial motion of H2 or ^He with respect to OCS. Sec- 
tion 5 summarizes and provides conclusions. 



2 Ground state 

The dimers 0CS-pH2 and OCS-^He are weakly bound van der Waals com- 
plexes which demand a full quantum mechanical calculation of the ^He degrees 
of freedom, due to their large zero-point motion. Exact methods [16,17,15] are 
available to determine ground and excited states of the dimers. However, we 
are interested in larger clusters 0CS-(H2)Ar and OCS-^Hejv, with special fo- 
cus on the influence of the superfluid properties of the solvation shell around 
OCS [18,19]. For this reason, in this paper we employ diffusion Monte Carlo 
(DMC) which is applicable to large systems of sizes on the order oi N = 100, 
and compare with the result of exact calculations for = 1 . 

The theory and implementation of DMC has been described extensively in the 
literature [20,21,22]. We employ here a mixed branching-weighting algorithm. 
Reliable interaction potentials V between OCS and "'He, and H2 (assumed 
spherical), respectively, have been calculated by Higgins et.al. [6,23]. Addi- 
tional interaction potentials are available for the OCS-He system [24,25]. Fig. 1 
shows cuts of the OCS-H2 and two OCS-He potentials along the minimum en- 
ergy path. In order to minimize computational cost in DMC, all potentials 
were expanded in Legendre polynomials and the resulting expansion coeffi- 
cients interpolated with splines. This is responsible for the small irregularities 
in the lower panel of Fig. 1, but we have verified that these do not affect our 
results. In this work we will not be concerned with the determination and 
quality of the potentials in regard to experimental measurements, and shall 
therefore start our investigation with the given Hamiltonian 

where mx — m(H2) and m(^He), B is the rotational constant of free OCS in 
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its ground state, and M its mass. d(px and d(py are the infinitesimal angles 
of rotation of the principal axis frame of OCS, r = |rx — ro| is the distance 
between atom and center of mass of OCS, and 9 the angle between the OCS 
axis and the location of the atom. We employ the rigid body diffusion Monte 
Carlo algorithm[26,27] in which the rotational degrees of freedom (p^ and (py 
are sampled by random walks in the angTilar variable taken in the principal 
axis frame of the molecule. Although it is not essential for a system of only 8 
coordinates (r^, rQ,ipx, ^y), of which only r and u = cos 9 enter in the potential 
calculation, we employ biased DMC here, samphng the state $o^r instead of 
$0) with the trial wave function 



^T{r, u) = exp air" + (ti + rii{u - Wi)^) log ■ 

+ {t2 + V2r\u-U2f (2) 

which is flexible enough also for > 1 complexes. The parameters of <I>t are 
chosen to approximately minimize the energy expectation value, and are given 
in table 1. This trial function has both angular and radial dependence, and 
thus it requires the full rotational importance sampled algorithm developed 
in Ref. [27]. Importance-sampled DMC yields exact values for the ground 
state energy and also for expectation values of operators commuting with 
the Hamiltonian. Unless otherwise specified, all expectation values of opera- 
tors not commuting with the Hamiltonian are computed here with descendent 
weighting [28] implemented according to the procedure for pure estimator of 
Ref. [29]. This corrects the mixed expectation values that result from impor- 
tance sampling {i.e., averages over the product $o$t) by a factor of $o/$r 
to obtain the exact expectation value over $q. This is particularly important 
for structural quantities. 



3 Methods for Rotational excitations 



Since DMC transforms the Schrodinger equation to a diffusion equation whose 
asymptotic solution is the ground state, DMC is often very efficient for the 
calculation of ground state properties while imposing only a small well con- 
trollable time step bias. The situation is quite different for excited states. 
Exact calculation of excitation energies is possible with the projection opera- 
tion imaginary time spectral evolution approach (POITSE) which is generally 
very time-consuming. [30,31] Alternatively, excited states can be calculated 
with fixed node approximations, but these generally introduce an unknown 
amount of bias by the choice of nodal surfaces[31]. The methods described 
below share the element of approximation with fixed node in the sense that 
they are approximative, but possess the advantage that it is very straightfor- 
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ward and can be implemented with the same efficiency as ground state energy 
calculations. 



3.1 Clamped coordinate quasiadiabatic approximation (ccQA) 

We summarize here the procedure which Quack and Suhm [12] introduced to 
calculate rotational levels of the van der Waals dimer (HF)2. The rigid-rotor 
Hamiltonian H^ot of the instantaneous configuration R at diffusion time t is 
obtained by diagonalizing the associated instantaneous incrtial tensor /. For 
OCS-Xtv, its elements in the body fixed frame of OCS arc 



with Ti = J2-y xj^l- Diagonalization yields principal values I a, Ib, and Iq, and 
hence the rotational constants A = h/8n^clA, B = H/Sti'^cIb, C = h/8TT^cIc- 
For N — 1 these should correspond to those of an asymmetric top for both 
He and H2 complexes, as seen experimentally [7,8]. For given total angular 
momentum J, the resulting rigid-rotor cluster Hamiltonian is 



This is diagonalized in a symmetric rotor basis, [32] resulting in eigenenergies 
ej^T, T — — J, . . . , J. Since ifrot depends parametrically on the configuration 
R, the energies ej^r are also parametrically dependent on R. As pointed out 
by Quack and Suhm in Ref. [12], they can therefore be regarded as providing 
effective centrifugal potentials, Vj^^(R) = ej^r, which may be added to the clus- 
ter Hamiltonian (1) to yield an effective Hamiltonian for different rotational 
states of the cluster as a whole: 



Because the Hamiltonians -f^(R) and Hj^t(R) differ only by the centrifugal 
potential terms, implementation of correlated sampling [33,34] on all systems 
Hj^T-(R) with a single random walk, but employing different individual weights 

is comparatively straightforward. Each system results in an excitation energy 
Ej^r, where these are the ground state energies of the respective effective 
Hamiltonians. The ground state energy Eq can be calculated inexpensively 
together with all excitation energies Ej^^- in the single DMC run. 



H,ot^AJl + BJl + CJl. 



(3) 



//j,,(R) = i/(R) + 16,,(R). 



(4) 
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3.2 Eckart frame ccQA (EccQA) 



The diagonahzation of the moment of inertia tensor / at every (imaginary 
time) instant t amounts to choosing the principal axis of / as the instanta- 
neous molecule fixed axis system, which thus moves with the motion of the 
generalized coordinate R. However, the axis system which separates rotational 
and small vibrational motion most effectively is the Eckart frame [13] and this 
is the preferred axis system for the derivation of rovibrational Hamiltonians. 
Therefore, we modify the above ccQA prescription in the following manner to 
accommodate the preferred Eckart reference frame. 

(1) Instead of diagonalizing the instantaneous inertial tensor I (which is of 
course independent of the axis system used) , we calculate the components 
lij in the Eckart frame and invert this matrix. 

(2) Instead of solving for the eigenvalues of the asymmetric top Hamiltonian 
(3), we solve for the eigenvalues ej^^ of 



This Hamiltonian is also used in the method described in the next section. 

For the determination of the instantaneous Eckart frame, a reference geometry 
needs to be specified. In the case of OCS-X, it can be specified by the position 
of X with respect to the OCS center of mass and molecular axis, Tref and 
cos^ref- As reference values of these coordinates, we choose the ground state 
expectation values r^f = (r) and cos^rcf = {cos 9). For OCS-H2 and OCS-'^He 
wc find these expectation values are given by Tref = 3.704A and 3.935A, and 
^ref = 105.7° and 108.4°, respectively 

3.3 Rotational constants from ground state averages (GSA) 

For semi-rigid molecules, rotational and vibrational degrees of freedom can 
be treated in a Born-Oppenheimer approximation in which the Hamiltonian 
is averaged over the rotational ground state. This is achieved by defining an 
effective rotational Hamiltonian i/f? 



frame axes given in Ref. [14], with interchange of y and z. Since for small 





(5) 
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amplitude vibrations I' differs only slightly from the pure inertial tensor /, we 
have followed Ref. [14] and also calculated (^o|(-^~^)i,j|*o)- For OCS-X, the 
effective rotational Hamiltonian (5) takes the form 

H^^, = AJI + + CJI + D{jjy + JyJ,) (6) 

where y is defined as the axis perpendicular to the OCS-X plane. The exphcit 
expressions for A, B, C, and D are given in Ref. [14]. Rotational energies are 
then obtained by expanding /Jj'r^ in the symmetric rotor basis and diagonal- 
izing the matrix. As noted already above in Section 2, we compute the exact 
ground state expectation values in Eq. (5) by descendant weighting of the 
importance sampling, in order to remove the trial wave function bias. 

By an ad hoc combination of ccQA and GSA, one can also calculate rotational 
constants from (approximate) excited state averages in a similar fashion and 
use these to go back and calculate rotational energies. The procedure is as 
follows. First, sample an approximate excited state corresponding to (</, r) 
defined by the effective Hamiltonian Eq. (4) of the ccQA approach explained in 
section 3.1. Second, calculate the rotational constants A, B, C, D as averages 
over these approximate excited states, i.e., j^j-\A\^ j^t) etc. This needs to 
be done also with descendent weighting to correct for the trial function bias. 
Third, diagonalize iJ^^ to obtain an updated estimate of the excitation energy 
Ej^T- We found that this ad hoc estimate turned out to be inferior to the ccQA 
results for both complexes, 0CS-(H2) and OCS-^He. 



4 Results 

We have applied the DMC-based approximate methods described above to 
the calculation of rotational excitation energies and rotational constants for 
0CS-(H2) and OCS-^He and compared with the exact energies which we have 
determined using the program BOUND [15]. To relate the exact reference 
energies to rotational constants, we fitted the BOUND rotational energies to 
the eigenenergies of the Watson A-reduced Hamiltonian. [36] In case of the 
GSA calculation, the resulting rotational parameters A, B, C, D are used to 
construct the effective Hamiltonian (6) which we then diagonalized to obtain 
rotational energies. 

In all calculations reported here, the imaginary time step value in the DMC 
walk was dt/h = 0.00011cm for 0CS-(H2) and dt/k = 0.00022cm for OCS- 
^He. 

Fig. 2 shows the respective ground state densities of H2 and "^He with respect to 
OCS, p(r, 9). These densities are obtained by simply binning the coordinates of 
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H2 and ^He in the body fixed frame of OCS in an unbiased DMC calculation, 
i.e. with trial function = Ij and then squaring the amplitude. The densities 
are qualitatively quite similar for the two complexes, with the helium density 
being somewhat more delocalized. In particular, we note that ^He has a non- 
neghgible probabihty to reside near the south pole of OCS (the oxygen end, 
at negative z). This leads us to expect that a semi-rigid approximation for 
OCS-^He might be problematic. 



4.1 Exact excitation energies from BOUND 

Fig. 3 summarizes the exact and quasi-adiabatic energy levels for OCS com- 
plexed with H2 (top panel), and for OCS complexed with "^He (middle and 
bottom panels). The exact energy levels are shown as horizontal bars. For the 
helium complex, we show the results from two different OCS-He potentials, 
namely the MP4 potential of Higgins and Klemperer [6] (middle panel) and the 
HHDSD potential of Gianturco and Paesani [24] (bottom panel). The ground 
state energy of OCS-H2 is -74.59 cm'^ and of OCS-^He is -16.38 cm'^ (MP4) / 
-15.85 cm~^ (HHDSD). The exact results for OCS-H2 show a very well-defined 
energetic separation between the low energy rotational states and the higher 
energy states associated with the vibrational motion of the H2 molecule in the 
OCS-H2 potential. For low values of J, the total angular momentum, these two 
different sets of states are separated by ~ 25 cm~^ This separation decreases 
as J increases and the rovibrational and purely rotational states are expected 
to mix, but even for J — 6, the highest value shown here, there is still a clear 
separation of ~ 5 cm~^ between the highest rotational sublevel and the first 
rovibrational level. 

These results indicate that OCS-H2 can be described as a rigid rotor and sug- 
gests that dynamic approximations based on such an assumption will provide 
a good starting point for e.g., analyzing energy levels in larger clusters. This 
situation changes considerably for the lighter OCS-He complex. With both the 
MP4 potential (middle) and even more so with the HHDSD potential (bot- 
tom) we find the energy separation between rotational and vibration levels 
greatly reduced. For small J, this separation is only ~ 8 cm~^ with the MP4 
potential, and further decreases to ~ 4 cm~^ with the HHDSD potential. This 
smaller separation between the purely rotational and rovibrational levels for 
the helium complex can be ascribed to the weaker interaction between OCS 
and ^He. 

The energetic spread of the rotational sublevels r is also different between the 
two systems, with the helium complex levels being compressed relative to the 
corresponding levels in the hydrogen complex levels: the lighter mass of H2 
results in a large rotational constant A compared to B and C, which leads 
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to a steeper increase of the rotational energy with increasing r = Ka — Kc, 
i.e. for large Ka, which approximately corresponds to the fast rotation of H2 
about a tumbling OCS. 

However, this smaller spread for the helium complex is clearly outweighed by 
the decrease in rovibrational energies. As a result of the smaller gap between 
rotational and rovibrational energies for OCS-He, a mixing of rotational and 
rovibrational levels is predicted at J > 4 for the HHDSD potential (bottom 
panel), and at J > 5 for the MP4 potential (middle panel). In comparison, 
for the OCS-H2 complex the mixing does not occur until J = 7. The onset of 
this mixing of rovibrational and rotational levels at small J values is another 
indicator of non-rigid, "floppy", behavior, seen here in the extreme for the 
helium complex. Therefore, from energy consideration, pure rotations cannot 
be distinguished from rovibrations at large values of the total angular quantum 
number J. 

Table 3 and 4 list the parameters of the Watson A-reduced Hamiltonian ob- 
tained by fitting the energy levels obtained from the BOUND calculations to 
the following Watson A-reduced form: 

Ha^I{B + C)3' + [A-^-{B + C)]3l +^-{B- C){Jl - J^) - Ajj' 

-5k[JI{JI - 3l) + {Jl - 3l)Jl]. (7) 

For comparison we show the corresponding experimental values of these pa- 
rameters as determined spectroscopically in Refs. [7] and [8]. These were ob- 
tained by fitting to energy levels reaching up to J = 6, but still lying below the 
rovibrational levels. These tables provide a reference for the rotational con- 
stants derived from the approximate methods, described below, and provide 
a measure of the accuracy of the three potentials [6,23,24] considered in this 
work. Both complexes are fit by asymmetric top parameters. 

4-2 Rotational energies from ccQA and Eckart modification (EccQA) 

Both the assumption of separation of vibrational and rotational motion that 
is central to the ccQA approach, and the resulting effective Hamiltonian (5) 
break down when the respective vibrational and rotational energies are of 
similar magnitude. Thus from Fig. 3 and the discussion of the exact ener- 
gies above, it appears that ccQA approach should break down at J ~ 4 — 5 
for the helium complex, and at J ~ 7 for the hydrogen complex. The GSA 
calculation is also expected to fail at these J values. Furthermore, neither of 
these approaches can yield the higher rovibrational states, where, in addition 
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to the rotation of the complex, the H2 or ^He atom also vibrates radially 
with respect to the OCS molecule. To obtain estimates for the energies of 
these excited states possessing both vibrational and rotational character, one 
would need to sample the various excited vibrational states for J = (e. 5'., by 
imposing nodal constraints on the trial wave function [37,12], and then add 
the centrifugal potential Vj,r(R) to the Hamiltonian (1) as described above 
to obtain vibrationally excited states with J > 0. Since this would introduce 
more approximations, we refrained here from trying to obtain the rovibrational 
excitation energies with ccQA beyond the pure rotational energies. 

The ccQA energies for OCS-H2 and OCS-He are shown as crosses in Fig. 3. 
They generally track the corresponding asymmetric top rotational energy lev- 
els, showing the greatest accuracy for the more rigid OCS-H2 complex, and 
larger deviations for the more floppy OCS-He complex. For comparison, we 
have plotted as a continuous dotted line the quadratic behavior of the en- 
ergy levels of free OCS, i.e., BJ{J + 1), with B = 0.2028 cm'^. [18] It is 
notable that while both complexes are asymmetric tops, the free molecule line 
does nevertheless track the lowest set of rotational levels for OCS-H2, while 
for OCS-He the lowest rotational levels lie considerably lower than the free 
molecule values. In Figs. 4a, 4b, and 5, we show the error of the ccQA and 
EccQA calculations of rotational levels with respect to the exact BOUND en- 
ergies, AEj,r = -Ej,r((E)ccQA) — (BOUND), for values of J where pure 
rotational levels can still be identified. Not surprisingly, the errors introduced 
by the (E)ccQA approximation are larger for the less rigid OCS-^He complex, 
see the larger y scale in Fig. 5, relative to Fig. 4a. Nevertheless, the deviations 
of (E)ccQA from the exact energies are still quite small in both cases, consider- 
ing that both dimers, especially OCS-'^He, experience large zero-point motion 
(see Fig. 2). Comparing now the original quasi- adiabatic approach (ccQA, 
closed symbols in Figs. 4a, 4b, and 5) with the Eckart modification (EccQA, 
open symbols), we see that the introduction of the Eckart reference frame 
does not improve the accuracy of the results. For completeness, we also have 
compared the relative errors for OCS-He when the two different potentials are 
used, MP4 and HHDSD (not shown). Since the HHDSD interaction is less 
anisotropic. Fig. 1, and allows for larger zero-point motion of the ^He atom, 
the (E)ccQA approximation yields poorer results with th HHDSD interaction 
than with the MP4 interaction, by up to a factor of two for J = 1 — 3. 

Fig. 6 shows the mean distance (r) of H2 and ^He, respectively, from OCS 
for the rotational excited states (J, r), obtained in the ccQA approximation. 
Although (r) spans a smaller range for OCS-H2 due to the larger rigidity, (r) 
rises more steeply with r = Ka — Kc, because of the high energies involved 
with the quantum number Ka (see the comment above about the larger spread 
in energies evident in Fig. 3). 
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4.3 GSA 



The four rotational parameters of the ground state effective Hamiltonian i/^^, 
eq. (6), are tabulated in Table 5. Comparison with the fits to the exact energy 
levels shown in Tables 3 and 4 show that using the pure inertial tensor / 
instead of the effective inertial tensor /' tends to improve the agreement with 
the values of the asymmetric top rotational constants A, B, C from fitting to 
the BOUND results. We note that, other than for the ccQA approximation, 
GSA rotational constants are sensitive to the choice of axis system, and the 
Eckart conditions may not be neglected. 

The deviations of the eigenenergies of H^^^ from the exact rotational energies 
are shown in Figs. 7a, 7b, and 8 for 0CS-(H2) and OCS-^He, respectively. 
Comparison with Figs. 4a and 4b shows that, for OCS-H2, the GSA devia- 
tions are of similar magnitude as those of the ccQA approximation, with GSA 
doing slightly better for smaller J than ccQA and vice versa for larger J. For 
OCS-^He, the deviations of GSA increase faster with J. Also, table 5 shows 
that, just like ccQA, GSA is less accurate for OCS-^He than for 0CS-(H2). 
There is a small, but finite probability for the ^He atom to be close to the oxy- 
gen side, along the OCS axis, i.e. OCS-'^Hc vibrating end-over-end. For these 
configurations, the rotational constants A, C{C'), and D diverge, leading to 
poor statistics in the DMC calculation. 



5 Conclusions 



We have compared the exact rotational excitation energies of the weakly 
bound Van der Waals complexes 0CS-paraH2 and OCS-^He with two approxi- 
mate methods, namely, the clamped coordinate quasiadiabatic approximation 
(ccQA) and ground state averages of the rotational constants (GSA), as well as 
with variants of each of these methods. The purpose of this survey is to assess 
which approximate methods will work best for larger complexes 0CS-(H2)jv 
and OCS-^Hciv, in the absence of non-stochastic exact techniques. 

While for the GSA approach, fulfilling the Eckart condition improves the re- 
sults considerably, we found no significant improvement of the ccQA results 
when the instantaneous inertial frame is replaced by the Eckart frame (Ec- 
cQA). This is good news for larger size clusters, since incorporating the Eckart 
condition brings with it the need to define a reference structure, which would 
become increasingly cumbersome and less tractable for growing N. We also 
found that incorporation of a centrifugal correction in the instantaneous in- 
ertial matrix in the GSA method (Section 3.3) actually tends to degrade the 
accuracy of the excitation energies. This is in agreement with the findings of 
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Ref. [14]. 



A significant piiysical feature tliat emerged from comparison of tiie complete 
energy spectrum of rotational and rovibrational states for the two complexes of 
OCS with H2 and with He, is the greater rigidity of the hydrogen complex, as 
evident from the larger separation between purely rotational and rovibrational 
states. This was evident from surveying the overall pattern of energy levels. 
Fig. 3. Experimental investigation of vibrationally excited states for OCS- 
H2 and OCS-^He would be useful, because their energies are most sensitive 
to the potential surface V{r, cos 9) at higher energies and larger r values. 
This differentiation in rigidity apparent from examination of the excited state 
spectrum is particularly interesting given the close similarity of the ground 
state densities for the two complexes (Fig. 2). 

A detailed projection for the breakdown of the assumption of separability of 
vibrational and rotational excitations in the multi-particle complexes OCS- 
Xjv, N > 1, is of course difficult, because exact results for all excited state 
energies are currently impossible to obtain for large A^. However, the projection 
operator imaginary time spectral evolution approach (POITSE) can provide 
exact energies given a good trial wave function. It yields the energies of selected 
rotational states, where the selection is determined by the choice of a projector 
acting on the ground state. POITSE calculations [11] of rotational excitations 
of OCS-^He^v for A^ < 20 show approximate agreement with ccQA results 
for A" < 5 and J — 1, with the ccQA results lying always somewhat lower 
than the POITSE results. This agreement with exact calculations indicates 
that the rigid coupling approximation underlying the ccQA approximation is 
good in this regime. For larger A^, the POITSE rotational energy saturates 
in accordance with experimental results, while the energies obtained in ccQA 
and GSA monotonically decrease for increasing A". This imphes that not all 
^He atoms in the first solvation shell can rigidly follow the OCS rotation, as 
discussed extensively in Ref. [19] and demonstrated explicitly in Ref. [11]. 
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Fig. 1. The interaction energy along the minimum energy path of the potential 
energy surface V{r,cos9) between H2 and OCS (top) [23] and ^He and OCS 
(bottom). For the latter we show both potentials used in our calculations, 
MP4 by Higgins et.al. [6] (squares are ab initio points from that reference) 
and HHDSD by Gianturco et.al. [24] (dashed hne). For usage in the DMC 
simulations, all potentials have been interpolated from the ab initio data points 
by splines in radial direction and by expansion in Legendre polynomials in 
angular direction, see Ref. [24]. 
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Table 1 

Parameters for the trial wave function (2) where r is in A. 





ai 


a 


h 


m 


Ul 


02 


t2 


m 


U2 




c 


0CS-(H2) 


-6.6 


0.79 


6.4 


-1.28 


-0.22 


0.529 


-1.0 


-0.71 


-0.22 


5.9 


1.59 


OCS-^He 


-5.9 


0.79 


7.68 


-0.32 


0.0 


0.529 


-1.0 


-0.57 


-0.26 


6.3 


1.78 
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Table 2 

Equilibrium configuration and root mean square deviations. 





(r)[A] 


(cos 9) 


Ar[A] 


A cos 6* 


0CS-(H2) 


3.704 


-0.270 


0.364 


0.158 


OCS-^He 


3.935 


-0.315 


0.478 


0.200 
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Fig. 2. Probability density p{r, 6) for H2 (top) and *^He (MP4 potential of 
Higgins et.al. [6], bottom) to be found at distance r from the OCS center 
of mass and at angle 9 in the OCS fixed coordinate system, obtained from 
unbiased DMC. The OCS molecule is oriented along the 2;-axis, with oxygen 
at negative z, and the x-axis is defined as any axis perpendicular to z. 
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Fig. 3. Top panel: excitation energies Eji of 0CS-(H2), obtained by BOUND 
and in the clamped coordinate quasi-adiabatic approximation; middle and 
bottom panels: Ej^^ for OCS-^He, using the MP4 potential of Higgins et.al. [6] 
(middle) and the HHDSD potential of [24] (bottom). For comparison, the 
dotted line shows the quadratic behavior of the free OCS rotational energy 
levels BJ{J + 1), with B = 0.2028cm-^ [18] 
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Fig. 4a. The error made in ccQA (filled triangles) and EccQA (open trian- 
gles), AEj^r = ^j,r((E)ccQA) - (BOUND), for the complex of OCS with 
hydrogen, '0CS-H2, J = 1, 2, 3. 
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Fig. 4b. Same as Fig. 4a for J = 4, 5, 6. 
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Fig. 5. The error made in ccQA (filled triangles) and EccQA (open triangles), 

AEj^^ = Ej,^((E)ccQA)-Ej,^(BOUND), for the complex of OCS with hchum, 
OCS'-^He, J = 1,2,3, using the MP4 potential of Higgins et.al. [6]. Note the 
larger vertical scale, relative to that for 0CS-(H2) in Fig. 4a. 



23 




Fig. 6. Mean distance (r) between H2 (top) and *^He (MP4 potential of Higgins 
et.al. [6], bottom), respectively, and the center of mass of OCS, for J = 0, . . . , 6 
and T = — J, . . . , J, in ccQA approximation. 
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Tabic 3 

Parameters of the Watson A-reduced Hamiltonian, eq. 7, for 0CS-(H2) by fitting 
to the energies obtained by BOUND, J = 1, . . . , 6, and to the experimental energies 
of Ref. [8]. 





MP4 


exp.[8] 


A 


0.7595 


0.7607 


B 


0.1997 


0.19996 


C 


0.1540 


0.15344 


Aj 


1.4 X 10-6 


2.4 X 10-6 


AjK 


1.7 X 10-^ 


1.56 X 10-"^ 


Ai^ 


3.1 X 10-^ 


10-3 


5j 


3.7 X lO-'^ 


7 X 10-^ 


5k 


1.23 X 10-^ 


1.4 X 10-^ 
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Tabic 4 

Parameters of the Watson A-reduced Hamiltonian, eq. 7, for OCS-*^He by fitting to 
the energies obtained by BOUND, J = 1, . . . , 4, for both, the MP4 and the HHDSD 
potential, and to the experimental energies of Ref. [7] . 





MP4 


HHDSD 


exp.[7] 


A 


0.4253 


0.4309 


0.44059 


B 


0.1833 


0.1946 


0.1835992 


C 


0.1198 


0.1191 


0.1221314 


Aj 


6.260 X 10-5 


2.244 X lO-'' 


3.1378 X 10-5 




-5.40 X 10-^ 


-8.40 X 10-5 


4.669 X 10-5 


Ak 


1.23 X 10-3 


2.50 X 10-3 


1.43663 X 10-3 


Sj 


1.89 X 10-5 


3.42 X 10-5 


1.12265 X 10-5 


Sk 


4.38 X 10-^ 


1.55 X 10-^ 


3.2667 X 10-^ 



26 



Tabic 5 

Parameters of the effective Hamiltonian ^1- {^)-> OCS-H2 and OCS-"^He 

(MP4 potential of Higgins et. al. [6] ) , 





0CS-(H2) [cm-i] 


OCS-^He [cm-i] 


A 


0.7520(2) 


0.42(3) 


B 


0.20318(5) 


0.206(8) 


C 


0.15556(7) 


0.12084(2) 


a 


0.15663(5) 


0.13355(3) 


D 


0.0016(6) 


-0.07(2) 
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Fig. 7a. The error made in the ground state average calculation, AEjj. — 

Ej,,(GSA) - Ej,, (BOUND), for the complex of OCS with hydrogen, OCS- 
(H2), J = 1, 2, 3. Both inertial matrices, / (filled circles) and /' (open circles), 
were used. 
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Fig. 7b. Same as Fig. 7a for J = 4, 5, 6. 
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Fig. 8. The error made in the ground state average calculation, AE'j^^ = 
Ej,,(GSA) - Ej,, (BOUND), for the complex of OCS with helium, OCS-'^He, 
J = 1,2,3 (MP4 potential of Higgins et.al. [6]). Both inertial matrices, / 
(filled) and /' (open), were used. 
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